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Abstract 

Stereoscopic spectral imaging is an observing technique that affords rapid ac- 
quisition of limited spectral information over an entire image plane simultaneously. 
Light from a telescope is dispersed into multiple spectral orders, which are imaged 
separately, and two or more of the dispersed images are combined using an analogy 
between the (x, y, A) spectral data space and conventional (x, y, z) three-space. Be- 
cause no photons are deliberately destroyed during image acquisition, the technique 
is much more photon-efficient in some observing regimes than existing techniques 
such as scanned-filtergraph or scanned-slit spectral imaging. Hybrid differential stere- 
oscopy, which uses a combination of conventional cross-correlation stereoscopy and 
linear approximation theory to extract the central wavelength of a spectral line, has 
been used to produce solar Stokes-V (line-of- sight) magnetograms in the 617.34 nm 
Fe I line, and more sophisticated inversion techniques are currently being used to de- 
rive Doppler and line separation data from EUV images of the solar corona collected 
in the neighboring lines of He-II and Si-XI at 30.4 nm. In this paper we develop an an- 
alytic a priori treatment of noise in the line shift signal derived from hybrid differential 
stereoscopy. We use the analysis to estimate the noise level and measurement preci- 
sion in a high resolution solar magnetograph based on stereoscopic spectral imaging, 
compare those estimates to a test observation made in 2003, and discuss implications 
for future instruments. 

Subject headings: instrumentation: spectrographs, methods: analytical, techniques: 
spectroscopic 



1. Introduction 

Spectral imaging in general, and solar spectral imaging in particular, suffer from a fundamen- 
tal problem in detector physics. Spectral images have three independent variables (x, y, A), while 
current image detectors only support two independent variables (x, y) and integrate over wave- 
length A. Conventional techniques to overcome this problem generally use time-multiplexing: in 
filtergraph imaging spectroscopy, a narrow band filter is tuned slowly across the spectral range 
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of interest and an image collected at each discrete A; in conventional scanned- slit imaging spec- 
troscopy, the light is passed through a spatial filter (the slit), selecting a single x, and the remaining 
light is dispersed to project A onto the detector's x axis. Even more sophisticated techniques such 
as Fourier imaging spectroscopy use time multiplexing to collect multiple two-dimensional basis 
images of the three-dimensional data space. Multiplexing in time is photon-inefficient as photons 
that are not part of the current sample are discarded. That is a problem because solar measurements 
are photon starved: instruments must race to collect sufficient photons for a spectral measurement, 
before the features on the Sun change. 

Stereoscopic spectral imaging overcomes the problem of spectral imaging by analogy be- 
tween the spectral imaging problem and the stereoscopic problem of determining feature distance 
in ordinary 3-space. Dispersed spectral images are integrals along diagonal lines in (x, y, A) space; 
they are analogous to images of a 3-space subject, with the images collected at a "look angle" 
that depends on the dispersion of the instrument. Collecting multiple dispersed images in dif- 
ferent spectral orders yields data that can be inverted stereoscopically to measure some spectral 
characteristics of the line everywhere in the image plane simultaneously. 



Kankelborg et al.l (|200ll) and iFox et al.l (120031) have described using this stereoscopic ap 



proach to simultaneous imaging and spectroscopy in EUV emission lines. Their MOSES rocket 
payload, launched 2006 Febr uary 8, obtained Doppler m easurements of explosive events and jets 
in the He II line at 30.4 nm (|Kankelborg and Foxl 120071) . Another related technique, computed- 
tomography imaging spectroscopy (CTIS) uses multiple spectral orders for hyperspe ctral imaging 



and spectropolarimetry at low spectral resolution i n the visible and infrared (e.g. [Wilson et al 



19971 : Miles 



1999; Miles et al.ll 19991 : lDereniakll2005l) . 



DeForest et al. (2004) have developed a first-order theory of stereoscopic inversion that is ap- 
plicable to a much simpler spectral context: measuring Doppler shift and the longitudinal Zeeman 
splitting in a narrow visible absorption line in the solar photosphere. That work is conceptually 
similar to the MOSES stereoscopic imaging, but there are two key differences: there is only a single 
spectral line in the instrument passband, and it is an absorption line; and the line width is narrower 
than a pixel in dispersion. DeForest et al. demonstrated the technique with proof-of-concept ob- 
servations of the quiet Sun, a decayed active region, and a clean sunspot. However, in that work 
they did not analyze the noise characteristics of the observing and inversion technique, only report 
the a posteriori noise measured in a quiet Sun region. 

In this article, we derive expressions for the a priori noise level (hence sensi tivity) and sys 



tematic error of a stereoscopic instrument in the narrow-line regime described by DeForest et al. 



(|2004|) using hybrid differential stereoscopy to determine line shift. We use a priori statistical 
calculation to derive the effect of the dominant source of uncorrelated noise - photon counting 
statistics - and also address quasi-random systematic errors in the inversion using some basic as- 
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sumptions about the scene being viewed. 

The analytic expressions are general and may be applied to stereoscopic spectrographs view- 
ing either absorption or emission lines, but we consider them in the specific context of an absorption- 
line Zeeman magnetograph viewing the relative line shift between the right- and left-circular po- 
larizations. In Owe analyze each of several noise sources individually; in §[3] the results of the 
individual analyses are applied to generate a noise budget for an example baseline instrument ob- 
serving the quiet photosphere with -0.1 arc sec resolution and 0.03 arc sec pixels,in the 617.34 nm 
Fe I absorption line. Finally, in §@] we draw conclusions about the types of observation for which 
hybrid differential stereoscopic spectral imaging offers the best prospects for advances over the 
state of the art. Throughout the discussion we have used A to indicate a single difference from the 
correct or expected value of the following quantity, and angle-brackets to denote RMS averages 
over an image. 



2. Hybrid differential stereoscopy 

Here we develop a theory o f noise in a hybrid differential analysis of absorption line data. 



Following iDeForest et al.l (|2004f) . consider an image in the (x, y) plane produced by observing, 
through a slitless dispersing spectrograph, a narrow absorption line of rest wavelength A , offset 
A'(x, y) from that wavelength, constant width AA, and integrated total intensity absorption L(x, y) 
over that width; limiting the wavelength range on the detector is a filter with transmission profile 
F(A') and total admitted continuum intensity C(x, y) in the absence of the spectral line. Without 
loss of generality, take the dispersion direction of the spectrograph to be along the x axis and 
consider the images collected in the ±1 spectral order. Then the y variation of the image drops out 
of the analysis and inversion may be performed independently along each line parallel to the x axis. 
Variations in the central wavelength of the spectral line result in spatial distortions of the images; 
these distortions are antisymmetric across the two spectral orders, and can be used to recover low 
spatial frequencies of the function A'(x) by correlation stereoscopy. Small patches of the images 
are cross-correlated to determine the offset functions X Lj ±i(x), and the wavelength offset can then 
be determined by 

2a 

where a is the dispersion of the +1 order. The L subscripts refer to low spatial frequencies; the 
maximum spatial frequency that may be resolved by this method is determined by the patch size 
used for the cross -correlation. The inversion used for Eq. [T]is analogous to the inversion carried out 
by the human visual system, determining the z coordinate of observed objects by cross -correlating 
the visual fields from the left and right eye. Recovering higher spatial frequencies requires a more 
sophisticated inversion. DeForest (2003) and DeForest et al. (2004) developed the process of 



-4- 



differential stereoscopy to determine high spatial frequencies in A' based on first-order intensity 
variation in the dispersed images. The difference between the +1 and -1 order dispersed images 
yields several terms in first approximation order, one of which is proportional to dA' /dx. By 
integrating the difference signal along the x direction it is possible to recover an approximation 
to the high spatial frequencies A' H throughout the image plane - because the primary detection 
term is to the spatial derivative, low spatial frequencies are not well measured with that technique. 
Furthermore, differential stereoscopy is sensitive to "leakage" of intensity and other signals into the 
dA' /dx measurement. These terms are proportional to A' itself and must therefore be compensated 
in some way in strong field regions. 

By resampling the dispersed images to eliminate the spatial shifts related to A' L before apply- 
ing differential stereoscopy, it is possible to eliminate image-related noise in the differential signal, 
to first approximation order. The differential process is used to recover high spatial frequencies, 
and the low spatial frequencies are inserted via a mathematical filter. The resulting symmetric 
hybrid inversion equation for first-order stereoscopic inversion (from DeForest et al. 2004) is: 



where A' L is again the low- spatial-frequency component of the central- wavelength shift image, as 
determined by explicit cross -correlation of patches of the +1 image with corresponding regions of 
the -1 image, X L> _i, and Xl,+i are again distorted spatial coordinates that use the correlation signal 
to remove measured offsets in the images, I ±1 is the intensity profile of each image as a function 
of focal-plane coordinate x or the distorted focal-plane coordinate X Lj ±i, a is the dispersion of 
the instrument, and e is a small convergence factor that splices the frequency spectra of the two 
derivation methods. 

Note that Eq. [2]is implicit rather than closed-form: A^is present both on the left hand side of 
the equality, and also inside the integral on the right hand side. The presence of the corrective filter 
term e(A' cl — A' L ) is what makes Eq. [2] a hybrid inversion: it preserves the high spatial frequencies 
present in the main portion of the integrand, while forcing the low spatial frequencies in the integral 
to match the A' L signal obtained by direct cross-correlation of image patches. If A' L happens to be 
constant and zero, then the impulse response function of the resulting spatial filter is a decaying 
exponential. In practice, the integral should be performed numerically in both the positive-X and 
negative-X direction, and the two results summed. The forward integral is causal along the x axis 
and the reverse integral is anti-causal, so that the resulting impulse response function is symmetric 
in x. 

The discrete form of Equation [2l then, is best calculated using the sum of two discrete filter 
transforms. Let J be the normalized difference image in compensated coordinates, interpolated to 
(integer) non-distorted pixel locations, and let i run over the pixel coordinates in the x direction. 
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Then J may be written: 

J z = 77 rVTl TvY (3) 

where the X subscript is to indicate that the I±\ images have been interpolated so that the i th 
pixel is found from the focal-plane location X\ ± aA L {xi). In the first-order theory developed by 
DeForest et al. 2004, J is important because, when dA'/dx <C a^ 1 , 

LdA ' 

which allows recovery of A' by integrating and scaling appropriately. 

The integral becomes a sum of two discrete spatial filters in the discrete symmetric hybrid 
inversion equation: 

^ - ^ (± ((!L ™ + <> - - 1 «!L J(n) - <{n) ) (1 

(5) 

where pixel centers are considered to exist at all integers between and w, md(C/L) RMS is the 
average ratio of admitted continuum intensity to line absorption intensity. The eA' L terms merely 
produce a smoothed version of A^: if A' L were constant, and the image were infinite in extent, then 
each of the two summations would return exactly A' L by the law of geometric sums. The low pass 
filter from this smoothing effect exactly matches the high-pass filter applied to the integral by the 
exponential decay term, so that the two versions of A' are combined together seamlessly. Equation 
|5j together with the process of cross-correlation to determine the relationship between Xl,±i and 
x, is a complete discrete image inversion. 

Figure[[]shows the step response of the filter described by Eq. [51 to illustrate how the different 
parts of the inversion work together. The step is located at position 100. The high-pass response to 
the differential stereoscopy signal (dark blue) is due to the interaction between the differentiation 
inherent in the detection method, and the "imperfect" (high-pass) integration imposed by the eA^ 
decay term in the integrand of Eq. |2l The exponential decay is matched by the low-pass response 
to the correlation stereoscopy signal (light cyan). The sum of the two signals is the value of the 
the RHS of Eq. |5j and reproduces the original signal. The slight overshoot of the reconstructed 
(magenta) curve is a discretization effect due to the phase of the Heaviside function relative to the 
pixel grid: a Heaviside function with the opposite phase relative to the pixel grid would yield an 
equal-magnitude undershoot. As with any discrete sampling filter, Eq. [5J operates best on spatial 
frequencies less than 2/3 of the pixel Nyquist frequency, because the reconstructed amplitude of 
those frequency components is independent of their phase relative to the pixel grid. The Heaviside 
function contains frequencies both below and above that limit, and the high frequencies' recon- 
structed amplitude depends on phase. 



-6 



In practice, errors in both the measurement of intensity and the inversion affect the analysis, 
so that the measured A' cl is not exactly the same as A' even if the image frequency spectrum is 
appropriately limited. In this error analysis we consider six principal sources of error: photon 
counting noise, (AI) RMS , that represents an additive, uncorrelated noise source at each pixel in 
each of the differential and correlation stereoscopy steps; misalignment noise (AXl) rms in the 
cross -correlation step; the effect of intensity gradients (both dC/dx and dL/dx) on the line shift 
signal; a nonzero slope dF/d\' of the instrument filter transmission function F(\') (the slope is 
assumed to be zero in Eq. |2]); and an excessively high dh! jdx signal in the solar image. The first 
two noise sources arise (directly or indirectly) from imprecision in the measurement of / in each 
pixel, and the last four sources arise from the non-ideal nature of the instrument and/or the solar 
field being imaged. In the following subsections we treat each error source independently. 



Photon counting noise is due to uncertainty in the measured intensity J ±1 in each focal plane, 
through the statistics of photon counting. Each measurement does not determine intensity directly, 
but rather fluence of energy (or, equivalently, photon count) on each pixel during an exposure time. 
Intensity can be derived by knowing the size of the pixel and the length of the exposure. Each pixel 
measures a fluence $ with an uncertainty AQ p h. given by: 



so that the RMS value (AI/I) RMS = (A$/$) HM5 = ^/hc/A^. This effect contributes to A' cl 
primarily via a photon noise term that must be propagated through Equation [5J the much higher 
photon counts in many-pixel image "patches" greatly reduce photon counting noise in the cross- 
correlation portion of the inversion ( §12.21) . The RMS photon noise contribution (A p hA' cl ) RMS is 
just the incoherent sum of the photon noise terms in each term on the RHS of Eq. [5J Each of those 
terms, in turn, is (to first order in A ph I / 1 and neglecting dl/dx for the calculation of noise in each 
pixel): 



where n p h is the per-exposure photon count in the corresponding pixel in each of the two orders. 
Eq. [7J simply reproduces the familiar behavior of incoherent counting statistics for intrinsic (nor- 
malized) quantities. 

The treatment may, at first glance, appear overly simple as the individual pixels in the images 
are interpolated and hence do not represent individual samples of the A ph & random variable but 



2.1. Photon noise in differential stereoscopy 
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rather scaled incoherent sums of adjacent samples. This is acceptable to ignore because, in the ab- 
sence of strong gradients in dX Lt ±i/dx (i.e. gradients comparable to unity per pixel), the samples' 
noise profile is preserved by the inversion. Adjacent pixels are co-added in the discrete inversion. 
Under linear interpolation, adjacent interpolated fractions of the signal from each detector pixel 
are approximately conserved by the interpolation and added coherently in the cross-pixel sum in 
Eq. [5J Because addition is commutative, the noise effect of the interpolated pixels on the final sum 
is the same as if no interpolation had taken place, up to a negligible factor of (1 — e). Put another 
way, although considering interpolation would further reduce the noise in Jj for individual pixels, 
co-adding Jj for adjacent pixels in the inversion eq. [5j would then include coherent cross-pixel 
noise terms that would exactly cancel the reduction - so we may safely choose, instead, to ignore 
the effects of the interpolation and use incoherent sums throughout. 

Carrying Eq. [7] forward through Eq. [5J yields: 

1 /C\ 
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Neglecting any field edge effects (in other words, taking < i C mj) permits taking the infinite 
sums directly, yielding (to first order in epsilon): 



Thus the statistics of the differential stereoscopy step are worsened by the width of the correction 
term. This is reflective of the fact that random walks diverge, and is the reason for using a hybrid 
inversion rather than a direct integral inversion: using direct integration in Equation [51 rather than 
high-pass integration with correction, is equivalent to setting e to zero, which would yield an 
arbitrarily large amount of photon-counting noise in the individual measurements. In actual use, it 
is necessary to balance e and the cross correlation patch size: larger patch sizes improve correlation 
behavior but require smaller e and thus worsen the photon noise due to integration of the differential 
signal. 

It is important to note that the photon noise from the discrete inversion is not independently 
sampled at each pixel: it is the sum of many noise terms along the x direction, so that the difference 
measurement between adjacent pixels has a much smaller noise level, yielding higher precision for 
small features. In particular, shifting by 1 pixel re-samples only one e th of the population of random 
samples in the calculated A' cl , reducing the photon noise by a factor of \fi so that small features 
are only subject to: 
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which is to be expected considering that the difference signal is the original signal being measured 
by the differential stereoscopy - the difference in measured wavelength between two adjacent pixels 
is simply a linearly scaled form of Jj. Differential stereoscopy is more sensitive to small features 
than to large ones, at least from the standpoint of photon statistics. 

Taking "typical" values of e = 0.09 (FWHM of 15 pixels), C/L = 10, and n ph = 5 x 10 4 
yields an absolute partial sensitivity a (A ph A' cl ) RMS = 0.052 pixel, and a gradient sensitivity 
a (A ph (A' cl) j — A'cij-i)) jims = 0-01^ ( me gradient sensitivity is unitless: pixels per pixel), due 
to photon counting alone. 



2.2. Photon noise in correlation stereoscopy 

The second source of error in Equation \5\ is misalignment in the cross correlation step that is 
used to derive X £j±1 from the original images. In the 2004 derivation of Eq. [2l it was necessary 
to assume that the displacement of the images was small, leading to dl/dx terms that enter in first 
order because the +1 and -1 order intensities are measured at slightly different spatial locations. 
Here, we consider the effect of small misalignments due to errors in the correlation co-alignment 
process. In the test measurements made by DeForest et al., the cross -correlation step used a full 
two-dimensional patch of image around each sampled point, rather than (as in the differential 
stereoscopic portion of the inversion) merely the particular horizontal line surrounding each point. 
This is reasonable in part because the image is, in fact, a two-dimensional field rather than a col- 
lection of independent rasters; and in part because the technique was found empirically to reduce 
correlation "misses" and improve the image signal. Here, we consider the potential sources of 
error injected by cross -correlation alignment of patches of image around each point. 

Because correlation co-alignment is a nonlinear process and depends on the field in the image, 
it is necessary to make some assumptions about the field being viewed. We here confine our 
analysis to magnetic measurements of the quiet Sun, which consists principally of nearly uniformly 
bright granules separated by dark, well defined lanes. Some other types of field (for example, 
the UV chromospheric network) are analogous and may be considered with the same or similar 
analysis, but other types of field (such as sunspot penumbrae) are likely too complex for simple 
analysis and require numerical modeling. 

In this analysis we consider cross-correlation of an idealized small patch of quiet Sun with 
TiTp pixels, of nearly uniform intrinsic brightness everywhere except in a fraction f d of the pixels of 
which are occupied by narrow, dark features with steep sides that are, on average, t d pixels across, 
and an absolute intensity deficit of I d (Figure . For typical quiet-Sun visible light observations, f d is 
about 0.1 and I d is about 0.5 C. Assuming that the correlation algorithm does not 'miss' (it achieves 
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Step response function: hybrid differential inversion 
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Fig. 1 . — Step response of the filter in Eq. [5] shows how the different components of hybrid 
inversion combine to reconstruct the wavelength shift signal. The source wavelength (black) is a 
discrete Heaviside function with step at pixel 100. The differential stereoscopy step (in the absence 
of noise) yields a high-pass filtered version of the Heaviside function (blue), while the correlation 
step yields a low-pass filtered version (cyan). Because the filters are constructed mathematically 
to match one another, the sum of the two filters (magenta) reproduces the step. The low-pass curve 
suffers from edge effects at the far sides of the image, but in a typical detector with more than 10 3 
pixel width, edge effects are minimal. 



Fig. 2. — Sample correlation field used for our correlation photon analysis. Correlation between 
two images is used to align two patches of pixels with radius r p . The image field consists of dark 
lanes against a background of approximately uniform brightness (left). The lanes consist of side- 
walls with thickness td and depth Id (right). A fraction fd of the image is considered to be lane 
wall. 
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a match between corresponding features in the opposite order images), we assume that errors in 
the resulting overall position of the patches is due to randomization errors from the displacements 
(real and apparent) of the dark features between the two images. In the worst non-pathological case, 
every intensity feature in the image patch is displaced by a different random amount, representing a 
strong high- spatial-frequency field. In that case, each feature's apparent location is the incoherent 
sum of two random variable samples: one due to photon counting statistics in the sloping sides 
of the feature, and one due to the random displacements of the features themselves. The former 
depends on the steepness of feature sides, the overall brightness of the image, and the number 
of pixels in each feature wall. The latter depends on the dispersion of the instrument and the 
variability of spectral offset across the image 

The uncertainty in location of each feature due to photon counting statistics is related to the 
slope at the sides of the feature: correlation fitting converts intensity in the sloping sides of a 
feature to position, so we can estimate the position/wavelength noise due to error in the intensity, 
by simply scaling the intensity noise to account for this conversion. The cross -correlation requires 
comparison between two features, so one can write: 



where Aids is the uncertainty in brightness in the side-of-feature pixels, and the right-hand ex- 
pression accounts for the typical number of pixels in the sides of a particular feature. The n p h 
comes from the fraction (AI) RMS // and is the number of photons counted during the exposure. 
The unitless 27 /(2/o — Id) arises because the average brightness in a feature sidewall is halfway 
between the background intensity and the dark feature intensity, The extra factor of 2 under the 
radical comes from the fact that the noise in offset is the incoherent sum of the location noise in 
each of the two features being compared. 

The ensemble average error position (AX ph ) RMS is then just the incoherent average of the 
offsets of each feature wall, or 



Taking the patch to be about 15 pixels across, and assuming a scene of solar granulation with f d = 
0.05, t d — 2, I d — 0.5I , and n ph = 5 x 10 4 yields a noise estimate of {AX ph ) RMS = 2.4 x 10 -3 
pixels. Photon noise in the correlation step is hence negligible compared to photon noise in the 
differential step. This is not surprising considering the large number of pixels that are averaged 
over by the correlation step. Put another way, the noise spectrum is dominated by high spatial 
frequencies, which are "smoothed out" by the correlation process. 
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2.3. Image noise in correlation stereoscopy 

There is an additional source of error in the correlation step: because the wavelength signal 
A'(x) includes high spatial frequencies (the signal that is being measured by the differential stere- 
oscopy step), different parts of the image must be shifted in different directions to find the average 
offset over the entire patch being correlated. These high spatial frequencies introduce an error into 
the measurement of the shifts X i ±1 (x), because the correlation step weights the importance of the 
different offsets according to the average sidewall slope in each portion of the patch being corre- 
lated, rather than uniformly throughout the patch of image. This error term is independent of the 
photon noise, and is introduced even for successful cross-correlations. It depends on the inequality 
of the weighting of the features, as well as the amount of displacement: 



5A' patch 5\dI/dx\ 

{AX wgt ) RMS = — — (13) 

where 5 is used to indicate actual differences from the patch- wide mean (different from A, which 
is being used to indicate statistical error in a quantity), A' patch is the signed difference between A' 
at each location in the patch being correlated and the mean for the whole patch, dl/dx is the slope 
of the intensity function at each location in the image, and S\dl/dx\ is the difference between the 
current value of the slope and its patch-wide average. As above, angle-brackets symbolize RMS 
averaging over a complete data set; here, an over bar is used to represent signed averaging over the 
patch. The inclusion of both 5-quantities (on the RHS) and A-quantities (on the LHS) is because 
Eq. [13] describes errors in the measurement of X that are due to the interpretation of the images 
themselves even in the absence of injected noise, rather than to an injected noise source in the 
direct measurement of brightness in each image plane. 

The numerator in the RHS of Eq. [13] depends only on the scene being viewed. The signed 

— 1/2 

average simply applies an n sa m P scaling of the overall error term, according to the number of 
independent samples n samp of image data in the patch. The quantity being averaged is the product 
of two random variables with zero mean and hence can be expanded simply in terms of their 
covariance C(x, y): 

(AX wgt ) RMS = " — (c (6A' patch , 6 \dl/dx\) + (8A' patch ) RMS (6 \dI/dx\) RMS 

\ al I ax l RMS V / V J 

(14) 

where n f is the number of "features" - separately resolved regions that each represents an inde- 
pendent sample of image data - in the patches being correlated. 

Assuming that the wavelength and brightness signals are uncorrelated eliminates the covari- 
ance term: 



/A y \ a ( 6 Katch) RMS (S\dI/dx\) RMS 

[^ wgt ,uncorr) RMS ~ (dI/d X ) RMS y/flf 
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so that, if wavelength and intensity of emissions be uncorrelated, the location uncertainty of a given 
cross-correlated patch is proportional to the RMS variation in central wavelength and the RMS 
variation in slope of the image intensity. For example, taking the patch to be about 15 pixels across, 
rif = n p /4 (one unique feature every four pixels), a (5A' patch } RMS = 0.1 pixel, (dI/dx) RMS =0.2, 
and (5 \dl /dx\) RMS =0A yields (AX wgtjUncorr ) RMS = 0.006 pixel. 

In practice, the image intensity and wavelength variation are not completely uncorrelated. For 
example, small magnetic features tend to be found in dark intergranular lanes, rather than randomly 
across the quiet Sun, so that C {5A' patchl 51) ^ 0. It is not clear, however, how or whether dl/dx 
and A' are correlated. In the worst case, they are proportional and the covariance term is equal to 
the product term in EqO doubling the effective noise level from this source to (in the example) 
0.012 pixel. In the best case, 5 A' would be near zero everywhere except where dl/dx is small. In 
that case, the covariance term would nearly exactly cancel the product term in Eq. Q3]and eliminate 
this noise source. 

In the context of a photospheric magnetograph, some further information about the covariance 
term may be used. In the quiet Sun for example, the strongest magnetic fields are observed to 
exist in the intergranular lanes, near local minima in continuum intensity. This suggests that the 
covariance between dl/dx and 5 A' tch is in fact negative in real scenes of solar granulation, so that 
Eq. [TBI likely overestimates the strength of this noise source in quiet sun magnetic measurements. 

In scenes that contain more contrast, particularly magnetic-correlated contrast, image corre- 
lation problems can dominate the Zeeman shift signal, in particular, near sunspots very strong 
variations in intensity are strongly correlated to very strong variations in magnetic field, reduc- 
ing the effectiveness of this technique for inverting stereoscopic data of sunspot fields. This was 
observed by DeForest et al. (2004), who found that correlation introduced strong artifacts into a 
not- well-resolved sunspot penumbra. 



2.4. Image brightness gradients 

Equation|2]includes several important engineered cancellations. Failure of those cancellations 
causes leakage of intensity or other information into the derived wavelength signal. The most 
important, and most likely to fail, cancellation is the elimination of an intensity spatial-derivative 
signal from the J introduced in Eq. [3j In the narrow-line approximation described by DeForest et 
al. (2004), the single-image brightness is given by (their equation 7): 

\E C ® F] (x/a n ) , . . x dA' , s dL ,dLdA' 

u „ Jj _ L(x) + anL(x) _ + M _ _ 2atK> (x) -- (16) 
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(where, again, E c is, the continuum brightness per unit wavelength, F is the prefilter function, L(x) 
is the total line depth intensity, and a n is the dispersion of the n th spectral order). The differential 
inversion arises from noticing that the difference between two opposite orders contains a useful 
dA' / dx term: 

I n (x) - J_ n (x) = 2a n L(x)^ + 2a n A'(x)^ (17) 

in which the dA' /dx term is desired and the dL/dx term is undesired. In hybrid stereoscopy the 
images are resampled into Xl coordinates with the substitution: 

X L ,n = x + a n A' L (x) (18) 

that eliminates the undesired dL / dx term in Eq. [F7]to first order. This process inevitably introduces 
some other undesired terms. In particular, substituting Eq. [T8]into Eq. [16] yields (writing Cl(x) 
for the convolution term): 



In(X L , n ) « (c L (x) + a n A' L ^j - (l(x) + anA' L ^J 
( r dA' n , , dL dA' p a / a / d 2 L 

,dL n . , . , d 2 L n . , dA' dL 



+ ( a n A'^ + a 2 n A'A' L — + a 2 n A' L 



dx " dx 2 n dx dx 
)lt dLdA' / / dA'\ 2 dL ll d 2 LdA! K ,dLd'A' 



n dx dx ^ n L y \ dx ) dx ^ dx 2 dx dx dx 2 

where C(x) is the convolution term from Eq. [161 As before, subtracting opposite orders removes 
all terms that contain an even power of a, and doubles the remaining terms. Grouping all terms 
and discarding terms beyond first order in A' yields: 

r /v n r /v \ n ( T dA' K ,dC L ., dL n2 fdA'\ 2 u dL\ 
UX L>n ) - J_ n (X L ,_ n ) « 2a n L— + A L —j— + A' H - - 2a 2 J — ) A' L - (20) 

where the first term is desired and the other terms are noise. The dC^/dx term is normally small 
both because the filter function is symmetric in wavelength and because the convolution integral 
smooths across several pixels. The A^ term represents leakage of the dL/dx signal into the data. 
(Recall, A' H is just the high spatial frequency component of the wavelength shift, given by A'— A^). 
The low-frequency image contribution to the noise has been reduced by the square of the (normally 
small) factor a n dA' / dx. 
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The LHS of Equation [20] is merely the numerator of the J definition in Equation [3] The 
denominator is the sum of the two intensities, which retains the even powers of a n while canceling 
the odd powers. Again, many cancellations occur, leaving: 

I n {X L>n ) + I„ n {X L>n ) « 2 (c L {x) -L + 2a 2 n A'A' L ^\ (21) 

so that the summed brightness of the resampled images has a systematic noise term that is pro- 
portional to the second derivative of the line depth. Combining the terms yields a more complete 
expression for J: 



+ A^ + (A i ,-(q„ ! f) 2 A ; 
(Cl(x) - L) + 2a5A'Ai£§ 



dL 
L I dx 



* + A^ + A- *L _ M^i^ [ 22) 



Cx I L Ldx \ H \ n dx J L ) Ldx Cj, dx 2 dx 

where first order approximations to the fraction have been taken, and squares of perturbations have 
been ignored. The first term of Eq. [22] reproduces the desired signal (Eq. [4]) and the additional 
terms represent noise sources due to the image itself. Writing them as an error term AJ grad , 

a n A' L ( dC L f A' ( dA'\ 2 \ dL t d? L \ 

where the first terms represent leakage of intensity information into the wavelength signal (both 
from continuum variation and from line variation), and the final term represents leakage of rapid 
brightness variations into the J gain and is negligible in nearly all cases. 

The continuum intensity leakage term is dependent on large-scale variations in continuum 
brightness, because C^is smoothed (by convolution with the filter function) over a distance of 
a n SX, where 5 X is the width of the entrance filter. In typical applications, a n 5\ might be 5- 
15 pixels. A scene of solar granulation with C/L-10, aA' L ~0.1 pixel and a n 5\~10 pixels, with 
intergranular lanes 3 pixels wide and continuum dips of 30% in the lanes, will have a continuum 
intensity leakage error of ~5 x 1CT 4 in J/I in the regions on either side of the lane, equivalent to the 
noise from 10 6 photons. The error accumulates coherently over half of the dispersed filter width, 
or about 5 pixels, yielding a peak injected noise signal in aA' of -0.025 pixel in the lane center 
from just the continuum term. 

The continuum intensity leakage is much stronger in regions with strong large scale variations 
in intensity, coupled with strong line shifts (such as found in sunspots). In a sunspot scene with 
C/L-10, aA' L ~ 4 pixels, and a continuum dip of 75% within the sunspot, pixels near the sunspot 
edge will have a continuum leakage term of 0.2 in J/I, much larger than any photon noise in the 
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system. The error accumulates coherently for the full dispersed width of the entrance filter, yielding 
integrated errors in aA' L as high as 20 pixels near the edge of the sunspot. This large error injection 
is the primary reason why DeForest et al. (2004) found sunspot measurements to be challenging 
with this inversion technique. 

The line intensity leakage term is largest in regions of very rapid variation of the central 
wavelength, but is negligible compared to the continuum leakage in most scenes. Considering a 
"challenging" scene of solar plage with aA' L ~ 1 pixel, aA' H ~ 0.3pixel amplitude, adA' /dx~ 0.1, 
C/L ~ 10, and dL/Ldx ~ 0.5/pixel yields intensity leakage of about 7 x 1CT 3 into J/I from the 
central term. This is equivalent to the photon noise from 3 x 10 4 photons, and if sustained across 
a full correlation patch of width 15 pixels could inject an error signal of up to -0.1 pixel into the 
final inversion; this represents about 20% of the background field. With a more typical quiet-Sun 
case of aA' L < O.lpixel (30 G), aA' H ~ O.lpixel, and adA' /dL ~ 0.02, the intensity leakage noise 
is under 10~ 3 in J/I and thus leads to errors of under 0.015 pixel after accumulating over a 15-pixel 
FWHM patch. 

The rapid brightness variation leakage depends on rapid changes in the line depth, in re- 
gions that are magnetized. It is negligible in nearly all cases other than sunspot penumbrae. 
In the quiet Sun, the strongest value of the term occurs where g-band bright points are present 
in intergranular lanes. In a typical g-band bright point, with aA' ~ lpix, aA' L ~ O.lpixand 
d 2 L/dx 2 ~ 0.005/ pix~ 2 , J can be expected to have a relative error of only 1 x 10~ 3 from this 
term. In sunspot penumbrae, where both aA' L and aA' may be larger than order unity, and where 
d 2 L/dx 2 may be a factor of 4-10 higher than in intergranular lanes due to the fine nature of penum- 
bral roll structure, J can be expected to have a relative error of order close to unity from this term, 
severely degrading results in sunspot penumbrae even if the continuum leakage signal (above) were 
not a problem. 

In all cases, the intensity gradient leakage terms are proportional to A' and therefore represent 
a gain error rather than an uncorrelated additive noise source that affects detection of lone magnetic 
features. 



2.5. Filter passband gradients 

Filter passband gradients enter because the overall measured intensity is affected by the 
throughput function of the prefilter. In conventional filtergraph instruments, the slope of the side- 
lobes of the main filter passband is used to convert wavelength offset into an intensity signal; in 
stereographic instruments, that signal is undesired and represents a source of noise. Even in a 
properly tuned stereoscopic instrument, large excursions of A' carry the spectral line out of the 
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flat central region of the prefilter's passband and into the sloping wings. This modifies Eq. [T6l to 
include the effect of the filter on the measured L(x): 

1Bt ? /a -' ~ L(x)F(A>) + a n L(x)F(A')^ ] 

+a n K\x){F(K')f x +m% d -t) (24) 
-2^A'(x)(F(A0f^ + L(x)f:(f) 2 ) J 

where F(Q) is taken to be unity. Of the two new terms, only one is asymmetric, so the difference 
equation becomes: 

/. - U = F(A' W ) (2a„L W (l + f + Wl'wg (25) 

where both the overall factor of F and the additional unitless derivative (A'/ F) (dF/ dA') represent 
perturbations on the gain of the measured dA' /dx, and the final term is the intensity leakage that 
is removed by the hybridization step in the overall inversion. In a well-tuned instrument, the 
maximum value of F occurs near A' = 0, so we may expand F and F' around that value: 

1 [A 1 - A' \ 2 
F « 1 - - f — (26) 



2 V W 

where A'j ilt is the difference between the center wavelength of the filter and the rest wavelength 
of the spectral line, and W is the half-width at half maximum of the filter passband. Then the 
derivative dF/dA' is given by (A« K — A')/W . As an example, if A' filt is about -W/10 and A'is 
about W/4, this yields a gain error of 7% in the dA' /dx signal. 



2.6. Rapid variation of central wavelength 

Rapid variation of the central wavelength of the observed spectral line affects the measured 
signal from differential stereoscopy. The differential inversion (Eq. [5]) is derived using a first-order 
expansion of the brightness in terms of the slope dA' /dx of the central wavelength versus position. 
If the condition dA' /dx <C a' 1 does not hold (the slope is significant compared to the reciprocal of 
the instrument's dispersion) then the differential inversion equation is not valid and the derived A^ 
will differ from the true A' even in the absence of any other noise source. This is due to the effect 
being exploited to extract the dA' /dx signal from the intensities: small shifts in the A direction 
project into the dispersed spatial direction, so that the line of integration in the (x, A) plane that 
contributes to a single pixel may have more or less contribution from the line core depending on the 
relative angle between a n and dA' / dx. The resultant variation in intensity is extremely nonlinear 
when a n and dA' / dx are nearly parallel. Further, the dA' / dx intensity variations are due to small 
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shifts in image position asA' varies. If A' grows larger than 1 pix/a n , then the image shift becomes 
noticeable and distorts the resultant line-shift images. 

The various effects of large excursions in A' on spatial distributions are illustrated in Figure 
|3l which illustrates the geometry of the (x, A) plane in a stereoscopic measurement of line center 
for a spectral line with structured Doppler signal. Three forms of perturbation on the line are 
shown: a sinusoidal modulation of the line center wavelength; a triangle wave modulation; and a 
dual step function. The sinusoidal modulation fits the criteria outlined by DeForest et al. (2004) 
for differential stereoscopic inversion: a n dA'/dx ^> 1 and a n A' is not large compared to the 
spatial features of interest. The others fail in different ways. The triangle wave fails because 
a n dA'/dx ~ 1, and the offset fails because a n A'is large compared to the width of the transition 
features at the edges of the step in central wavelength. 

Hybrid differential stereoscopy reduces the effects by eliminating the low spatial frequency 
components, so that the amplitude and slope of the high spatial frequency component A' H , rather 
than of A' itself, are of import: structures larger than roughly the width of the patch used for 
correlation are attenuated or eliminated by the high-pass filtering that generates A' H . 

By considering carefully the geometry of the differential inversion in the (x, X) plane it is 
possible to estimate the systematic errors due to slope and line-displacement effects alone on dif- 
ferential stereoscopy, even in the absence of image effects such as variations in line or continuum 
intensity. Figure 0] shows the geometry of dispersed integration through a spectral line. The contri- 
bution L' n to the measured intensity I n is approximated to first order in Eq. |H but the contribution 
may also be written exactly in the case of constant L(x), constant dA' /dx, and a broad filter: 



I n (x) = J dX (F(X')E c (x - a n X') - L5(X - A'(x) + (a n X) dA'/dx) = C-- 



L 



+ a n (dA'/dx) 
(27) 

Applying this expression to Eq. [3]gives an analytic formula for the relation between J and dA'/dx. 
Letting 7 = a n dA'/dx and £ = L/C, 

J " = xrf = " 7 (rn) G-wa-o) ) <28) 

where the factor a n £/(l — £) is a linear calibration coefficient on dA'/dx and the right-hand term 
is a nonlinearity term. The nonlinearity in J n becomes strongly evident when 7 grows to about 
0.5, where (with L/C = 0.1) the nonlinearity term has a value of 1.38, and grows to infinity where 

7=v / 0^)- 

The strong nonlinearity in J is partially compensated by spatial distortion. The nonlinear term 
causes growth in J compared to the proportionality with dA'/dx, but that growth is accompanied 
by spatial distortions that partially cancel the nonlinearity. In particular, the measurable quantity of 



-18- 



interest is A' rather than dh! / dx. In a pure differential inversion, with 7 positive and non-negligible, 
the J n signal arises from a slight weakening of the spectral line in the — n spectral order and a strong 
strengthening in the +n order. However, the image in the +n channel is strongly foreshortened 
while the image in the — n channel is weakly fore-lengthened: the total ^-integrated intensity in 
each spectral order remains the same under small to moderate (<1 pixel) A' perturbations in the 
image. 

Consider a finite shift in central wavelength A' as illustrated in Figure @J Accompanying 
the shift in brightness is a spatial distortion that affects the two spectral orders differently. In the 
positive order the line depth is more strongly expressed, because it is concentrated into fewer pixels 
than might be expected from the spatial extent of the feature. Likewise, in the negative order the 
line depth is more weakly expressed, because it is spread across more pixels than the usual spatial 
extent of the feature. The corresponding variation of J n thus has three regimes in this simple 
system (neglecting the line width): 

{I a; I > 5x +^ 5A 
5x-a n SA I I ^ c _ a n 8A j 
2 di h -> <jx 2 j mi d — Io+I _ nl (zy; 
5x— a n SA ^ | | J _ In,i—I-n,i 

2 ^ l X l Jctr ~ I n ,l+I-n,l 

where (from Eq. [27]) I Q = C - L, I_ n j = C - L/(l + 7), and I njl = C - L/(l - 7). Simple 
analysis shows that 

— 7 

J ™ d = (i + 7 )(2/£-l)_l (30) 



and 



Jctr ~ (1 - 7 2 )/£- 1' (31) 



so that the measured wavelength shift may be computed directly: 



C f 

SA meas , = -"n 1 ^ / Jndx (32) 
= -a-H- 1 haJA (- ) + Sx(l - 7) ( ' ' 



(2-7)(l-l)- T ; v -V(l- 7 2 )-^ 
27^A \ (5xa?)(l->y)(>y) \ 

(2- 1 ){l-l)- 1 ) l-£- 7 2 J 
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where the last row expands the 7-dependent error factor to first order in 7. Equation [34] shows that 
the linear approximations used to derive [2] are quite good, at least for this simple geometry: the 
first nontrivial term in the error factor is two orders of approximation smaller than the extracted 
signal (taking both 7 and t to be small, as did DeForest et al. 2004). The high quality of the 
approximation can be seen in Figure [3] the measured peak and feature-averaged AA meas in each 
feature remains close to the injected AA of the model, despite gross violation of the assumptions 
used to derive Equation [2] 

The astute reader will have noticed that Equation[33]has a pole at 7 = \f\ — 1. That is because 
we have neglected finite length effects in the region where the line shift occurs. Equation[33]is more 
complex but exact for this geometry, provided that the line width is small compared to 5x/a n . 
When that condition is violated (which happens when the center locus is exceptionally narrow on 
the detector), individual pixel lines of integration do not integrate across the entire foreshortened 
spectral line but are limited instead by the ends of the spectral feature This limitation keeps the 
signal response J ctr finite. Equation [33] neglects this effect, leading to a spurious pole. 

For nonzero line widths, Equation [331 reveals a constant correction factor of 1/(1 — £) that 
must be applied to stereoscopic inversions in general. This correction (of order 10%) is a true 
calibration coefficient, in the sense that it is independent of 7. 



3. Summary & Specific Observing Scenario 

The effects in §[2] above separate into two primary classes: additive and multiplicative noise 
in the data. The additive effects, such as photon noise, affect sensitivity for detection of small 
magnetic or Doppler features; the multiplicative effects, affect the precision of magnetic measure- 
ments. Here we calculate a noise budget for a plausible observing scenario using a stereoscopic 
instrument as a photospheric magnetograph. 



Consider an instrument observing quiet sun thr ough the Dunn S olar Telescope (DST; iDunn 



1969 ) and the medium-order adaptive optics system (Rimmele 2004), at the Fe I 617.34 nm line, 
with 0.07 arcsec pixels (the diffraction limit is 0.2 arcsec), a 67 pm prefilter passband width detuned 
from A by 6.7 pm (10% of the bandwidth), £ = 0.1, dispersion a = 100 pixel ran -1 , n p h = 5 x 10 4 
photons per pixel, and a correlation patch FWHM of 15 pixels (e = 0.09). Consider that a "typical" 
resolved feature strength is between lOOG-lkG, with length scales of a few pixels (so that only 
stereoscopic measurements apply, without the stabilizing influence of the more precise correlation 
stereoscopy). The 617.34 nm line has a Lande g factor of 2.5, so that the Zeeman splitting equation 
is just ' 

he 

B' = ——(S'A) = b.GTnm^S'A = 5.6 kG ^S'A (35) 
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Fig. 3. — Five plots showing arbitrary variations of spectral line central wavelength show the 
manner in which the approximation of Equation \5\ breaks down at high line slopes. In order: the 
X-Lambda plane; the X-lambda plane, dispersed in two opposite directions (Pixel brightnesses are 
formed by vertical integration in these two plots); the original and derived spatial derivative of 
line center wavelength; and an overlay of the original and reconstructed signals. The sinusoidal 
segment at left is well matched. The pathological sawtooth and square wave segments are distorted 
by the coupling between wavelength and spatial position. 
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Fig. 4. — A finite linear feature in A demonstrates the interplay between projection, brightness, 
and feature size in a stereoscopic instrument. At top, a spectral line undergoes a finite shift SA 
between the points x' and x[ in the (x', A) plane. The instrument collects two spectral orders with 
dispersion a± n . Each order integrates pixels along lines of the appropriate slope: positive (cyan) 
or negative (red) and projects onto the detector x coordinate. At bottom the line depth brightness 
integral — (/„ — J_ n ) varies as a function of pixel location (projected back to the A = A line at 
top). 
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where B' is measured resolution- smoothed magnetic field strength ("flux density"), h is Planck's 
constant, c is the speed of light, g is the Lande factor, fj, is the Bohr magneton, and 5' A is the 
splitting of the line (twice the shift). The splitting measurement is the result of two independent 
measurements of a shift AA: one in right-circularly polarized light (Stokes I+V) and one in left- 
circularly polarized light (Stokes I-V). Hence the noise AB' is a factor of v2 higher than might be 
supposed by simple scaling of the A A noise terms in §[2] with Equation [351 

The noise budget for a quiet-Sun DST observation may be summarized as in -Table [Q Two 
principal types of error enter into the measurement: random or quasi-random noise as described in 
§ §12.1H2.4l and calibration nonlinearities that enter as a result of filter detuning and the first-order 
stereoscopic inversion as described in S S2.5H2.6l 

Because of the multiple sources of error, calculation of error budgets for actual observations 
must be carried out on a case-by-case basis; however, we note that for magnetograph observations 
of the photosphere, the dominant noise term is (as expected) photon noise for a plausible instrument 
setup and a range of targets that includes most non-sunspot-related targets on the Sun. The expected 
formal single-exposure magnetogram sensitivity is -10 Gauss (twice the RMS background noise 
level) and the expected magnetic flux calibration accuracy is expected to be close to 10% in most 
cases including plage. 

The error budget in Tabled] is broadly consistent with test measurements collected by DeFor- 
est et al. (2004) at the Dunn Solar Telescope. In those measurements, exposures of only ~8x 10 3 
photons were collected; 16 exposures were combined to create each image for an effective expo- 
sure of 1.2 x 10 5 photons; and an asymmetric inversion (whose sensitivity was reduced by a factor 
of \/2 compared to the symmetric inversion discussed here) was used. The resulting calculated 
photon noise level is thus a factor of a/24/5 higher than shown here, or 9.5 G (0.45 G for the 
correlation noise), and the non-photon noise sources are scaled up a factor of y/2 compared to the 
sample measurement described here. This leads to a total a priori noise level estimate of 10 G for 
that measurement. Indeed, in a small quiet region with no detected flux features, DeForest et al. 
found a posteriori RMS noise levels between 9-12 Gauss, in good agreement with the current a 
priori calculation. 

4. Conclusions 

We have calculated a detailed noise budget for hybrid differential inversion of the line-shift 
signal from a stereoscopic spectrograph, with specific application as a stereoscopic magnetograph 
observing a specific absorption line (Fe I 617.34 nm) in the solar photosphere. The noise budget is 
photon dominated and careful a priori analysis shows that it offers similar calibration performance 
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to existing and p ast filtergraph (SOHO/MDI, Scherrer et al.ll 1995 ) and Fourier tachometer (GONG; 
Leibacher 1999 ) instruments ( Jones and Ceja2001 ): instrument response to field is within a few 



to ten percent of linear over a wide range of pixel- averaged field strengths. 

The principal difference between the stereoscopic technique and the filtergraph or Fourier 
tachometer techniques is that the instrument can operate much more rapidly, because all required 
photons are collected simultaneously in a dual-polarized-beam instrument. This is very important 
for ground based observations in which atmospheric seeing effects severely limit the precision with 
which line shifts can be measured in time-multiplexing instruments such as filtergraphs. 



Some instruments, such as ZJJVIPOL (IPovellll998|) . have been able to limit crosstalk of atmo- 
spheric effects by multiplexing extremely rapidly; but such instruments are still subject to motion 
blur even in the absence of severe crosstalk. A stereoscopic spectral imager should be able to ac- 
quire a diffraction limited magnetogram with few-Gauss sensitivity in a small fraction of a second 
at most major observing facilities. We expect this rapidity of acquisition to be come important even 
in space-based observing as higher resolution instruments (e.g. Hinode/SOT: \lchimoto et al.ll2006l) 
release data to the scientific community: the race between solar evolution and data acquisition 
(photon counting) give rise to a tradeoff between sensitivity and spatial resolution for any spectral 
imager, even above the atmosphere. 

Particularly for features that are spatially small, stereoscopic imagers offer a surprising ad- 
vantage: because the fundamental signal extracted from the images is proportional to the spatial 
derivative of the line shift signal, such instruments are particularly sensitive to small features, 
right down to th e spatial resolu tion of the measurement. Like Dolby™} noise reduction for audio 
recordings (e.g. IVaseghill2006l) . spatial derivative measurement enhances high spatial frequencies, 
boosting the high frequency signal relative to the photon noise floor. Conventional spectral mea- 
surements sample each location in the image plane independently, so that in a roughly circular 
feature of linear size L and area A the uncertainty due to photon counting in the feature's total 
magnetic flux scales as A~ l l 2 or L^ 1 . In a differential stereoscopic measurement, the average 
slope uncertainty, rather than the total flux uncertainty, scales as L^ 1 so the total flux uncertainty 
is independent of feature size for features smaller than the correlation patches described in Section 
® 

Stereoscopi c spectral imaging is a special case of the more general technique of spectral to- 
mography (e.g. IWilson et al.lll997|) . and appears to offer significant advantages in the regime we 
have described: observations of a narrow spectral line, where only one or two moments of the 
line are desired to be measured. Spectral tomogr aphy using regularized inverse or conjugate anal- 
ysis techniques (e.g. IClaerbout and Fomell 12004) could potentially be used to extract much more 
information about the spectrum, using multiple orders or, as did Wilson et al., crossed gratings; 
but as the number of additional channels grows, so do the potential paths for signal contamination 
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and crosstalk between different aspects of the measured spectrum. We conclude that stereoscopic 
techniques represent a sweet spot for imaged Doppler or magnetic measurements of solar scenes, 
offering superior photon efficiency and comparable performance to filtergraph measurements. This 
is already important as ground-based spectral imagers are already photon- starved near the dif frac- 
tion limit in existing meter- class telescopes (such as the Dunn Solar Telescope, |Punnl ll969. and 
the Swedish Vacuum Tower. IScharmer et al.ll2003|) and will grow more so near the diff raction limit 



of fu ture telescopes such as the four-meter Advanced Technology Solar Telescope (IKeller et al. 
2002h . 
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Source 


Eq. /Section 


AARMS (pm) 


Afi'RMS (G) 


cal. err. 


Photon noise (diff.) 




0.53 


4.2 




Photon noise (corr.) 




0.024 


0.2 




Image noise (corr.) 




0.12 


1.0 




Continuum leakage 




0.25 


2.0 




Line leakage 




0.15 


1.2 




Filter signal leakage 


12573231 






3%-7% 


Stereo nonlinearity 








1%-10% 


TOTAL (in quadrature) 




0.62 


4.9 


3%-12% 



Table 1 : Noise budget for a sample stereoscopic magnetograph measurement as described in the 
text. This measurement would have ~ 6 G RMS noise even in field-free regions of the Sun, and 
3%-12% calibration errors in the measured flux of small magnetized regions of 100-1000 G, no 
more than a few pixels across. 



